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Abstract. In this paper we present an extension of standard iterative 
splitting schemes to multiple splitting schemes for solving higher order 
differential equations. 

We are motivated by dynamical systems, which occur in dynamics of the 
electrons in the plasma using a simplified Boltzmann equation. Oscilla- 
tion problems in spectroscopy problems using wave-equations. 
The motivation arose to simulate active plasma resonance spectroscopy 
which is used for plasma diagnostic techniques, see [2J, [TS] and [18] . 
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1 Introduction 

We motivate our studying on simulating a active plasma resonance spectroscopy 
which is well established in plasma diagnostic techniques. 

To study the model with simulation models, we concentrate on an abstract 
kinetic model, which described the dynamics of electrons in the plasma by using 
a Boltzmann equation. The Boltzmann equation is coupled with the electric field 
and we obtain coupled partial differential equations. 

The paper is outlined as follows. 

In section[2]we present our mathematical model and a possible reduced model 
for the further approximations. 

The functional analytical setting with the higher order differential equations 
are discussed in section [3] 

The splitting schemes are presented in in Section S) 

Numerical experiments are done in Section [5] In the contents, that are given 
in Section [51 we summarize our results. 
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2 Mathematical Model 

In the following a model is presented due to the motivation in [2], [16] and [15] . 

The models consider a fluid dynamical approach of the natural ability of 
plasmas to resonate in the near of the electron plasma frequency u) pe . 

Here we specialize to an abstract kinetic model to describe the dynamics of 
the electrons in the plasma, that allows to do the resonation. 

The Boltzmann equation for the electron particles are given as 

gfe^ = -v ■ V x f(x, v, t) - — V x cf> ■ V v f(x, v, t) 
at ra e 



-<j(x,v,t)f(x,v,i) + K,(x,v,v')f(x,v',t) dv', (1) 
Jv 

f(x,v,Q) = f (x,v), (2) 



and boundary conditions are postulated at the boundaries of P (plasma). 

In front of the materials we assume complete reflection of the electrons due 
to the sheath /(un + vj_) with vu is the parallel and v± perpendicular to the 
surface normal vector. d> is the electric field. 



3 Higher order differential equations 

We consider the abstract homogeneous Cauchy problem in a Banach space X € 
IR": 

, d n u(t) , d™- x u(i) , , , 

^ w _! =u«-i > t = l 1 ...,n, (4) 

where Aq, . . . , A n € X x X are bounded operators and || • || is the corresponding 
norm in X and let || • ||l(x) be the induced operator norm. 
For the transformation we have the following assumptions: 

Assumption 31 1.) The function fit) is given as: 

f (t) = 0, (5) 

otherwise we solve a non- autonomous equation. 

2.) We assume that the characteristic polynomial: 

A X n + A^- 1 + . . . + A n = 0, (6) 

has solution of complex valued matrices in X x X G C m x C m , given as: 

(XI - B 1 )(A/ - B 2 ) + . . . + (XI - B n ) - 0, (7) 
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Corollary 1. The higher order differential equation f3|) can be decoupled with 
the assumptions \31\ to the following differential equation: 



dui(t) 

dt 
du 2 (t) 

dt 

du n (t) 
dt 



- B lUl = 0, 

- B 2 u 2 = 0, 

- B n u n = 0, 



(8) 

(9) 
(10) 

(11) 



where the analytical solution is given as: 

n n 

u(t) = ^2 d i u i-o = X! ex P(Bit) diiiifi. (12) 

i=l i=l 

and di are given via the initial conditions. 

Proof. The solutions can be derived via the characteristics polynomial (idea of 
scalar linear differential equations) and the idea of the superposition of the linear 
combined solutions. 

Remark 1. The initial conditions are computed by solving the Vandermode ma- 
trix, see the ideas in [T7] . 
We have to solve: 



/ / 


/ . 


■ I \ 


Bx 


B 2 ■ 


■ B n 




Bl ■ 


■ Bl 


Ur 1 


sr 1 ■ 


■Bl- 1 ) 



( diui, \ 

^3^3,0 
\d n U nfi j 



( "(0) \ 
ffi«(0) 
&«(0) 



(13) 



A further simplification can be done to rewrite the integral-differential equa- 
tion in two first order differential equations. Later such a reduction allows us to 
apply fast iterative splitting methods. 

Corollary 2. The higher order differential equation $3$) can be transformed with 
the assumptions \31\ to two first order differential equation: 



dm(t) 



= B lUl (t) 



dt 

ui(0) = c?iui ;0 , 



du n (t) 



= B n u n (t) 



dt 

u„(0) = d„u„ j0 , 



(14) 

(15) 
(16) 

(17) 
(18) 
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where we have Bi — Bn + B^ for i = 1 . . . , n. 
The analytical solution are given as: 

n 

u(t) = ^2exp(Bit)diU i: o. 

»=i 

Proof. The analytical solution of the first order differential equation (|f4|) and 
(If 7p are given by each characteristic polynomial: 

Ai/-(Bn+Bi a )=Q 1 (19) 

(20) 

A„J - (B nl + B n2 ) = 0, (21) 

while the solution is given as with the notations: 

\iI = Ba+B i3) i = l ) ...,n, (22) 

and therefore the analytical solution is given as (fT9|) . 

Therefore this is the solution of our integro-differential equation ([3]) with the 
assumptions [3TJ 



4 Splitting schemes 

The operator-splitting methods are used to solve complex models in the geophys- 
ical and environmental physics, they are developed and applied in |22j . [23] and 
|24j . This ideas based in this article are solving simpler equations with respect 
to receive higher order discretization methods for the remain equations. For this 
aim we use the operator-splitting method and decouple the equation as follows 
described. 

In the following we concentrate on the iterative-splitting method. 



4.1 Iterative splitting method for Integro-differential equations 

The following algorithm is based on the iteration with fixed splitting discretiza- 
tion step-size t, namely, on the time interval [t n ,t n+1 ] we solve the following 
sub-problems consecutively for i = 0, 2, . . . 2m. (Cf. [15] and [H].) 



= Bacij{t) + B i2Cij ^{t), with cy(t n ) = d iC n (23) 
and Cio{t n ) = c n , c 4 ,_i = 0.0, 

dCl ^ t l(t) = B ilCij (t) + B i2 c iJ+1 (t), (24) 

with Cij+i(t") = diC n , 

where c" is the known split approximation at the time level t = t n . The split ap- 
proximation at the time-level t = t n+1 is defined as c n+1 = X)fc=i c k,2m+i{t n+1 )- 
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Theorem 1. Let us consider the abstract Cauchy problem in a Banach space X 
d t c{t) = B ll c{t) + B l2 c{t), 0<t<T 

i=l,...,n, (25) 

c(0) = diCi fi 

where Bn, B^, Bn + Bi2 :X — > X are given linear operators being generators of 
the Co-semi-group and c^o G X is a given element. Then the iteration process 
i.31\) I Sty) is convergent and the and the rate of the convergence is of second 
order. 

Proof. The proof is done in the work of Geiser P3] . 

The algorithm is given as: 
Algorithm 41 

= Ac l {t) + Bc l -\t), with c l (t n ) = c 1 - 1 ^ 1 ) (26) 
and the starting values c°(t n ) = c(t n ) results of last iteration , c _1 (t") = 0.0, 

dt {1 =Ac i (t) + Bc t+1 (t), (27) 
with c t+1 (t n ) = c\t n+1 ) , 

e > \c t+1 (t n+1 ) - c i-1 (i B+1 )|«op criterion (28) 

result for the next time-step (29) 

c(t n+1 ) = c m (t n+1 ), for m fulfill the stop- criterion (30) 

for each i — 0,2, ... , where c n is the known split approximation at the previous 
time level. 

In the following we concentrate on the iterative-splitting method. 

4.2 Iterative splitting method for higher order differential equations 

The following algorithm is based on the iteration with fixed splitting discretiza- 
tion step-size t, namely, on the time interval [t n ,t n+1 ] we solve the following 
sub-problems consecutively for j = 0, 2, . . . 2m. (Cf. [15] and [TT].) 

= B ilCij (t) + B i2Cij ^(t), with cy(t") = d lC n (31) 
and a (t n ) = c" , c 4 ,_! = 0.0, 

dCl ^ {t) = B a cij(t) + B t2 c lJ+1 (t), (32) 
with aj + i(t n ) = did 1 , 

where i = 1, . . . ,1 are the number of equations. Further c n is the known split 
approximation at the time level t = t n . The split approximation at the time- level 
t = t n+1 is defined as c" +1 = £Li c fc , 2m +i(* n+1 ). 
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Theorem 2. Let us consider the abstract Cauchy problem in a Banach space 
X C C: 

d t c(t) = B n c(t) + B i2 c(t), 0<t<T 

i = l,...,n, (33) 

c(0) = diCo 

where a\ G C is the constant based on the initial conditions, further Bn, Bi2, Bn + 
Bi2 : X — > X are given linear operators being generators of the Co-semi-group 
and Co G X is a given element. Then the iteration process \31 \) ~ nM) is convergent 
and the and the rate of the convergence is of second order. 

Proof. The proof is done in the work of Geiser [14] . 

The algorithm is given as: 

Algorithm 42 

Br (t) 

— ^ = B lti C i:j (t) + ll,,r, l With d(t n ) = &-\t n+1 ) (34) 

and the starting values Ci^o{t n ) — Ci(t n ) results of last iteration , Ci._i(t") = 0.0, 

dCi ^ t l{t) = B u c id {t) + B 2ti c iij+1 (t), (35) 
with c t+1 (t n ) =c l {t n+l ) , 

e > \c l+1 {t n+1 ) - c , '- 1 (* n+1 )|Stop criterion (36) 

result for the next time-step (37) 

c{t n+1 ) = c m {t n+1 ), for m fulfill the stop- criterion (38) 

for each j = 0, 2, . . . , where c n is the known split approximation at the previous 
time level. 

Further Bi = B\^ + B 2 .i is a decomposition of the matrix Bi. 

We reformulate to an algorithm that deals only with real numbers and 
rewrite: 

d t (c re (t) + ic lm (t)) = (Bre^i + iB im .n)(c re (t) + ic im (t) (39) 

+ (B re ,i2 + iBi m ,i2){Cre{t) + (t)), o < t < T, 

(c re (0) + ic im (0)) = (di tre + id iMl ){c refi + ic imfi )i = 1, . . . , n, 

We have the following algorithm: 

Algorithm 43 In the following, we have two iteration processes: 

— First iteration process: j — 1, . . . , J iterates over the decomposition of the 
matrices B\,B 2 . 

— Second iteration process: k = 1, . . . ,K iterates over the real and imaginary 
parts. 
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We start with j = 0, k = 0, 
First we iterate over j 

^(°i,Ve)(*) D 7".fe /,x . t-> J-l,k/.\ 
^ =-Bl,re<VeW + B 2,re< re : (*) 

- Bl,im<i,im 1 ( t ) ~ B 2 ,im,C i i ,im k (*)) 

— = Bl,re<' im (*) + #2,reC- )im ' (t) 

+ Si )im c|~e' fc_1 (*) + B 2,im4re' fc_1 (*)« 

witft ffte initial condition c{'r e (t n ) = cj'^ , c J i 'i m (t n ) = cj'^ 
with the starting condition c~^ k (t n ) — 0,c~^(t n ) = 



?"+i,fe/,\ 

= S 1)re d' re (<) + #2,re< re (*) 



0f 



— = B hre cil im {t) + B a ,r«<im (*) 

witt £/ie imft'ai condition 4+* ,fc (i n ) = 4tl' k ' n '> 41"™* = 41™*'" 



if j = J or the iteration error over j is less err then we iterate over k. 
Further c" is tte known split approximation at the time level t = t n , cf. Q/. 
Further Bi — + is a decomposition of the matrix B%. 

5 Experiments for the Plasma resonance spectroscopy 

In the following, we present different examples. 

5.1 First Example: Matrix problem with integral term 

We deal with a simpler integro-differential equations: 

dr r 

— - Ac(t)dt = B c{t') dt', t e [0, 1], (40) 

and the transformed second order differential equation is given as: 

d tt c =Ad t c + Bc (41) 

(42) 
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and the operators for the splitting scheme are given as: 



i—$6- 

while A T is the transposed matrix of A. 
The matrices are given as 

/-O.Ol 0.01 
0.01 -0.01 
0.01 0.01 -0.02 ... 
0.01 0.01 0.01 -0.03 



A A T 



- B 



(43) 



A = 



0.01 0.01 0.01 0.01 
0.01 0.01 0.01 0.01 



0.01 0.01 0.01 0.01 
V 0.01 0.01 0.01 0.01 

/-0.08 0.01 0.01 0.01 0.01 0.01 
-0.08 0.01 0.01 0.01 0.01 0.01 



B 



0.01 
0.01 



-0.08 


0.01 
0.01 



(44) 





-0.08 J 

0.01 \ 
0.01 



(45) 



-0.02 0.01 0.01 
00000 -0.01 0.01 
\ 00000 0.01 -0.01/ 

The Figure [T] present the numerical errors between the exact and the numer- 
ical solution. Here we obtain results for one-side and two-side iterative schemes 
on operators A and B. 

The computational results are given in the Figure [2j we present the one-side 
and two-side iterative results. 

The Figure [3] present the numerical errors between the exact and the nu- 
merical solution for the optimized iterative schemes. Here we obtain results for 
one-side and two-side iterative schemes on operators A and B. 

Remark 2. For the computations, we see the benefit of the optimal iterative 
schemes, which applied the two iterative steps of the two solutions in one scheme, 
see Section HI The best results are given by the one-side iterative scheme with 
respect to the operator B. 



5.2 Second Example: Third order differential equations 

We deal with a simple third order differential equations: 



< L£-Ac{t) = o, te [o,i], 
c(0) = (i,...,i)*eC"\ 



c'(Q) 



A x / 3 c(0), 



="(0) = \a^c(Q), 



(46) 
(47) 

(48) 
(49) 
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A e C m x C"\ c : IR+ C m is sufficient smooth (c e C 3 (1R + )) and we have 
m = 10. 

The transformed first order differential equations are given as: 



d t d - A 1/3 c x = 
d t c 2 -A^(-^ +l ^)c 3 =0 

d t c a -A 1 ' 3 {-^-i^)c 3 =0 



(50) 
(51) 

(52) 
(53) 

where c = diC{(t) and di, . . . , d 3 are given with respect to the initial con- 

ditions and are given as di = d 2 = d 3 = |c(0). 

Further the operators for the splitting scheme for the three iterative splitting 
schemes are given as: 



Ai^re = diag(A 1/3 ),Ai t 2, re outerdiag(A 1/3 ), 
M,i,re = -^diag{A 1/3 ),A 2 .2,re = -^outerdiag(A 1/3 ), 

M,\,im = —diag(A 1/3 ), A 2i 2,im = —outerdiag(A 1/3 ), 



A 3 ,l,re = -—diag(A 1/3 ),A 3a ,re 



The matrix A is given as 



^diag{A 1 ' 3 ),A 3Xi 



—outerdiag(A 1 ^ 3 ), 

■ —outerdiag(A 1 ^ 3 ), 



(54) 
(55) 

(56) 

(57) 

(58) 



(59) 



/-0.01 0.01 
0.01 -0.01 
0.01 0.01 -0.02 ... 
0.01 0.01 0.01 -0.03 ... 

0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 -0.08 
\ 0.01 0.01 0.01 0.01 0.010.010.010.01 -0.08/ 

Here, we deal with the following splitting schemes: 

ci is computed by a scalar iterative scheme. 

C2 7 c 3 are computed by a vectorial iterative scheme (because of real and 
imaginary parts). 

For c\ we have: 

dC\(t) 



dt 



= AuC[(t) + A 12 C\- L (t), with C\(t n ) = Ci" 1 ^ 1 ) (60) 



and the starting values C { {(t n ) = -C(t n ) 



(61) 
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where Ci = {c\. re + ici ( j TO ) and 



■*u=(V"„° )^ 12 =(V e /„ ), (62) 



Ai t i t re J 1 \ ^1,2,, 

for i = 1, 2, . . . , J and the solution is given as C\(t n+1 ). 
For c 2 we have: 

f)C i (t) 



dt 



= A2iC\{t) + A 2 2C\-\t), with C\{t n ) = C 2 ~ l (t n+1 ) (63) 



and the starting values C° 2 {t n ) = -C(t n ) (64) 
where C 2 = (c 2 , re + ic 2 , im )* and 

( ^2,1, re ) _4 — j ^2, 2, re ~ (^-2, l,im + ^-2,2 : im) Y,^ 

\ A2,l,re/ ' \ (^2,1, mi + ^2,2, im) ^2,2, re / 

for £ = 1, 2, . . . , J and the solution is given as C 2 (£" +1 ). 
For C3 we have: 

f)C i (t) 



Of 



= A 31 Cl(t) + A 3 2C l 3 - L (t), with CUt n ) = Cl- l (t n+l ) (66) 



and the starting values C° 3 (t n ) = -C(t n ) (67) 
where C 3 = (c 3jre + ic 3jml )* and 

_4 31 — [ ^3,1, re \ ^ 2 _ j A% : 2,re — (^-3, Mm + ^-3,2 : im) Jg^\ 

\ ^3,1, re / ' \ (^43 : 1,oti + ^3, 2, mi) ^3,2, re / 

for i = 1, 2, . . . , J and the solution is given as C 3 (t" +1 ). 
The solution is given as: 

c , (f" +1 ) = E 3 l 1 q-(f +1 ). 

The computational results for the optimized iterative schemes are given in 
the Figure 21 we present the one-side and two-side iterative results. 

Remark 3. For the computations, we see the benefit of the optimal iterative 
schemes. While we deal with real and imaginary parts, it is important to reduce 
the computational costs. We applied in one scheme the real and imaginary so- 
lution, see Section |4j The best results are given by the one-side iterative scheme 
with respect to the operator B. 



6 Conclusions and Discussions 

We present the coupled model for a transport model for deposition species in a 
plasma environment. We assume the flow field is computed by the plasma model 
and the transport of the deposition species with a transport-reaction model. 
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Based on the physical effects, we deal with higher order differential equations 
(scattering parts, reaction parts, etc.). We validate a novel splitting schemes, 
that embedded the real and imaginary parts of the solutions. Standard iterative 
splitting schemes can be extended to such complex iterative splitting schemes. 
First computations help to understand the important modeling of the plasma 
environment in a CVD reactor with scattering and higher order time-derivative 
parts. In future, we work on a general theory of embedding the complex schemes 
to standard splitting schemes. 
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Fig. 1. Numerical errors of the one-side Splitting scheme with A (upper figure), the 
one-side Splitting scheme with B (middle figure) and the two-side iterative schemes 
with 1, . . . , 6 iterative steps (lower figure). 
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Fig. 2. The computational time ol the one-side and two-side Splitting scheme: one-side 
splitting over A (upper figure), one-side splitting over B (middle figure) and two-side 
splitting scheme alternating between A and B (lower figure) with 1,...,6 iterative 
steps. 




Fig. 3. Numerical errors of the one-side Splitting scheme with A (upper figure), the 
one-side Splitting scheme with B (middle figure) and the two-side iterative schemes 
with 1, . . . , 6 iterative steps (lower figure). 
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Fig. 4. The computational time ol the one-side and two-side Splitting scheme: one-side 
splitting over A (upper figure), one-side splitting over B (middle figure) and two-side 
splitting scheme alternating between A and B (lower figure) with 1,...,6 iterative 
steps. 



